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Abstract 

A mostly right-handed sneutrino as the lightest supersymmetric particle (LSP) is an 
interesting dark matter candidate, leading to LHC signatures which can be quite distinct 
from those of the conventional neutralino LSP. Using SModelS vl.0.1 for testing the 
model against the limits published by ATLAS and CMS in the context of so-called Sim¬ 
plified Model Spectra (SMS), we investigate to what extent the supersymmetry searches 
at Run 1 of the LHC constrain the sneutrino-LSP scenario. Moreover, we discuss the most 
relevant topologies for which no SMS results are provided by the experimental collabora¬ 
tions but which would allow to put more stringent constraints on sneutrino LSPs. These 
include, for instance, the mono-lepton signature which should be particularly interesting 
to consider at Run 2 of the LHC. 
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1 Introduction 


Before the start of data taking at the LHC, the common perception was that supersymmetry 
(SUSY), if it has anything to do with stabilizing the electroweak scale, would be discovered 
quickly, while Higgs physics would need to wait for rather high statistics. In reality, quite 
the opposite has happened: a Higgs boson has been found (^[^, but there is still no sign of 
SUSY—or of any new physics beyond the Standard Model (BSM) whatsoever. 

Indeed, the searches at Run 1 of the LHC at centre-of-mass energies of 7-8 TeV have pushed 
the mass limits of SUSY particles quite high already, well above 1 TeV for lst/2nd generation 
squarks and gluinos [^|^. Scenarios with high-scale [^, split [6]-[^ or at least spread SUSY 
are thus becoming increasingly popular in the literature. It should be kept in mind, however, 
that the current LHC limits sensitively depend on the presence of particular decay modes, 
and are considerably weakened in case of compressed or stealth spectra. Besides, the 
squark/gluino mass limits vanish completely in case the neutralino LSP is heavier than about 
600 GeV. 

It should also be kept in mind that the SUSY mass limits depend sensitively on the nature 
of the LSP. Most experimental analyses indeed assume that the LSP is the lightest neutralino of 
the Minimal Supersymmetric Standard Model (MSSM). A particularly interesting alternative, 
and the subject of this paper, is a mainly right-handed (RH) mixed sneutrino in the MSSM 
augmented by a RH neutrino superfield 12,13 . This case is well motivated by two basic 


problems: the origin of neutrino masses and the nature of dark matter (DM). Its LHC signatures 
can be quite distinct from those of the conventional neutralino LSP. 

The left-handed (LH) sneutrino of the MSSM is excluded as the LSP and as a DM candidate 
because it has a non-zero hypercharge: its couplings to the Z boson makes it annihilate too 
efficiently in the early Universe, and hence its final relic abundance is lower than the value 
measured by the WMAP and Planck satellites [^[^. Very stringent limits come 
moreover from direct DM detection experiments: the vj^ scattering off nuclei is mediated by 
t-channel Z boson exchange, giving a spin-independent (SI) cross section of order 10“^®cm^ — 
a value excluded already a decade ago for DM particles heavier than 10 GeV. A light with 
mass below mzl2 is also excluded by the Z invisible width. 

The picture changes dramatically if we include in the MSSM a RH neutrino superfield 
(MSSM-I-RN from here on), which gives rise to Dirac neutrino masses. Besides the RH neutrino, 
the superfield also contains a scalar field, the RH sneutrino N (strictly speaking this is a 
right-chiral field, but we use the RH notation for simplicity). This field, if at TeV scale, can 
mix with the LH partner vl and yield a mostly RH sneutrino LSP as a viable thermal DM 

The phenomenology of this model was investigated in detail in in Indirect detection 


candidate 


and cosmology were discussed in (^1 [^[^, and LHC signatures in |25| ( |27| (see also |28| - [^ 


for related LHC studies). Reference M in fact gave an update of the status of the sneutrino 
as DM after the Higgs mass measurements, by exploring the SUSY parameter space with the 
soft breaking terms fixed at the grand unification (GUT) scale, and assessing also the impact 
of the most recent exclusion bound for DM direct searches from LUX IsTl. 


^Pure right-handed sterile sneutrinos can also be viable (non-thermal, depending on the model) DM candi¬ 
dates, as discussed e.g. in 16 - 20 . 
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In this paper, we extend the work of 27 by investigating to what extent the results from 


SUSY searches at Run 1 of the LHC, published in terms of so-called Simplihed Model Spectra 
(SMS) limits]^ constrain the sneutrino-LSP scenario. Moreover, we discuss the most promising 
topologies for which no SMS results exist but would enhance the LHC sensitivity to sneutrino 
DM. To this aim, we make use of the SModelS vl.0.1 package to compare the pre¬ 

dictions of the MSSM-hRN model against the SMS limits published by ATLAS and CMS. The 
strengths of SModelS vl.0.1 are that it 1.) automatically decomposes the signal of an arbi¬ 
trary SUSY spectrum into all its SMS-equivalent topologies, and 2.) includes a large database 
of more than 60 SMS results from ATLAS and CMS SUSY searches. This allows us to test 
the limits from a large variety of searches and at the same time draw conclusions about which 
additional topologies should be considered. 

The paper is organised as follows. After briefly deflning the MSSM-I-RN in Section we 
describe the numerical procedure in Section In particular in 3T we explain the sampling 
method and the constraints implemented in the model likelihood function, while in |3.2| we 
describe the application of SModelS vl.0.1 to the MSSM-I-RN. Our numerical results are pre¬ 
sented in Section]^ and the conclusions in SectionTwo appendices contain some interesting 
supplementary material. Appendix discusses the validity of applying SMS results from slep- 
ton searches (dilepton signature) to chargino-pair production followed by decays into leptons 
and sneutrinos. Appendix gives some details on scenarios with long-lived heavy charged 
particles, in particular gluinos or stops, which so far cannot be constrained by SMS results. 


2 The MSSM+RN model 


We use the MSSM-I-RN model as defined in (^,13,^. (The model used in 21,22,26 differs 
only slightly in notation.) The superpotential for Dirac RH neutrino superfleld is given by 


M- = e„(pHrH'‘ 




rIJ-, 


( 1 ) 


where Yj'^ is a matrix in flavor space (which we choose to be real and diagonal), from which 
the mass of neutrinos are obtained as Note that lepton-number violating terms 

are absent in this scheme. The additional scalar fields contribute with new terms in the soft- 
breaking potential 

U,oft = Lf Z/ + (M^)"^ N^*N^ - Hfl^NY + h.c.], (2) 

where both matrices and 1^/ are real and diagonal, Mjj = diag(m^fc) and = diag(Ap), 
with k = e, being the flavor index. In the sneutrino interaction basis, defined by the vector 

= 

(z>2, A*), the sneutrino mass potential is 


T/fc = _ *T 
^ mass 2 LR 




LR ^LR , 


with the squared-mass matrix Ai 


LR 


■^LR — 


-|- cos(2/5) -|- mZ 

^ ApU sin P — jiniD /tan P 


sin P — jirriD /tan P 


m 




+ mPjj 


( 3 ) 


( 4 ) 


^Simplified Models are effective-Lagrangian descriptions involving only a small number of new particles. 


They were designed as a useful tool for the characterization of new physics, see e.g. 32 33 . 
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Here, are the soft mass terms for the three SU(2) leptonic doublets, tan/3 = v^/vd and 
-\- v\ = (246 GeV)^, with Vu^ the usual Higgs vacuum expectation values (vevs). The 
Dirac neutrino mass m^) is small and can be safely neglected. 

The off-diagonal term determines the mixing of the LH and RH helds. If that 

is if the trilinear term is aligned to the neutrino Yukawa, this term is certainly very small as 
compared to the diagonal entries and is therefore negligible. However, can in general be a 
free parameter and may naturally be of the order of the other entries of the matrix 12,13 , thus 


inducing a sizable mixing among the interaction eigenstates. The sneutrino mass eigenstates 
are then given by 




— sin 6p cos 
cos 9^ sin 9r, 


n 

Nk 


( 5 ) 


The relevant parameters at the electroweak (EW) scale for the sneutrino sector are the two 
mass eigenvalues and and the mixing angle 6*^, related to the Ap term via 


sin 291 = V2 


Ap V sin /3 
(m? — m? ) 


( 6 ) 


The sneutrino coupling to the Z boson, which does not couple to SU(2)/, singlets, is largely 
reduced by a sizeable mixing. This has a relevant impact on the sneutrino phenomenology, as 
discussed in, e.g.^ Refs. (I^[^[^[W[[^ . 

The renormalization group equations (RGEs) are modified by the new singlet superhelds N 
as 


drra^fc 

dln/i 

dm|fc 

dln/i 

iAt 

dln/i 

dm^ 

_ 

dln/i 


(aY , 


IGtt^ 

(MSSM terms) 
2/3 


( 7 ) 


IGtt^ 


(Aif. 


IGtt^ 

(MSSM terms) 


-Agl _ Igl + 3 v 2 ^ ly 2 
2 ^ 10 1 2 * 2 ^ 


Ak 


E 

k=e,iJ,,T 


kA 


IGtt^ 


(Al) 


with /i being the renormalization scale, g 2 and gi the SU(2) and U(l) gauge couplings, Yt^j. the 
top and r Yukawa respectively. Notice that the RH soft mass receives corrections only from 


the trilinear term, which affects as well the running of the LH part, as recognized in 21,26|. 


By neglecting all lepton Yukawas but Yr in the RGEs and by assuming common scalar 
masses and trilinear couplings for all flavors, the sneutrino tau, ends up to be the lightest 
one among the three sneutrino flavors and hence the LSP, while z>ei = 9^^^. Note that it 
frequently happens that the mass splitting between and is smaller than 5 GeV, which 
means that regarding collider phenomenology they are practically degenerate. This will be 
discussed in more detail in Section 13.21 
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Table 1: Summary of the observables and constraints used in this analysis. 



Observable 

Value / Gonstraint 

Ref. 

Measurements 

rrih 

125.85 ± 0.4 GeV (exp) ±4 GeV (theo) 

1,^ 

) 


BR{B X,7) X 10^ 

3.55 ± 0.24 ± 0.09 (exp) 


39 



BR(R, ^ fi+fi-) X 10® 

3.2l};2 (stat) (sys) 


40 




0.1186 ±0.0031 (exp) ±20% (theo) 


15 


Limits 

y\pinvisible 

Z 

< 2 MeV ( 95% CL) 


41 



BR(h ^ invisible) 

< 20% (95% CL) 


42 




> 85 GeV (95% CL) 


43 




> 101 GeV (95% CL) 


41 



m-g 

> 308 GeV (95% CL) 


44 



^SI 

^ n 

^ ^Lux (90% CL) 


31 



3 Numerical procedure 


3.1 Sampling method over the model parameters 

For definiteness, we study the MSSM+RN with soft terms defined at a high scale M ~ 

Allowing for non-universalities in the gaugino and scalar sectors, our set of free 


27 


as m 
parameters is 

Ml, M2, M3, rriL, tur, mpf, uiq, mn, Ai,Ac,, Ag, tan/ 3 , sgn^u . (8) 

Here the M* are the gaugino masses, mL,mR, rriN are the charged slepton and sneutrino masses 
(equal for all flavors), mg is a common squark mass parameter, mn = 'i^Hu = denotes 
the common entry for the two Higgs doublet masses, and Ai and Ag are the scalar trilinear 
couplings for the sleptons and squarks respectively, same for all flavors. The absolute value of 
fi is obtained from the minimization of the Higgs potential, leaving only the sign of /r as a free 
parameter. The computation of the mass spectrum follows that explained in 27 , where all 
details are provided. 

The list of constraints implemented in the model likelihood function is given in Table In 
particular, besides consistency with H-physics constraints, we require the Higgs mass ruh to be 
compatible with the ATLAS and CMS measurements [^[^, which we combine by a statistical 
mean, as obtained in 45 . Its uncertainty is dominated by the theoretical error, estimated to 


be around 4 GeV 46 . We also require that chargino and charged slepton masses fulfill the LEP 


bounds at 95% confidence level (CL) —notice that the tau slepton has a slightly less stringent 
lower bound of 85 GeV (43] as compared to selectrons and smuons— and we include the gluino 


mass bound from the DO collaboration |44|. If is light enough to be produced in Z decay, we 


require its contribution to the Z invisible decay width to be smaller than 2 MeV 47 . Similarly, 


when the sneutrino mass is lighter than m/j/2, the Higgs can decay invisibly into sneutrino 
pairs. We require that such decays do not contribute more than 20% to the Higgs invisible 
branching ratio (4^. 
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Regarding DM constraints, we require consistency with the measured relic abundance and 
with the bounds from direct detection experiments (constraints from indirect DM detection are 
also fulhlled). The experimental error on has become incredibly small due to the Planck 

measurement , while the theoretical one is still large. We use a conservative estimate of the 
order 20% 48 for the latter. Furthermore, we enforce the sneutrino SI scattering cross section 


off nuclei, , to be compatible with the recent 90% CL bound from LUX 31 


To evaluate the experimental observables we first computed the supersymmetric particle 
spectrum with a modified version of SoftSusy [^. For the computation of the sneutrino relic 
density and elastic scattering cross-section the model has been implemented in FeynRules 


50 


51 , by adding the appropriate term in the superpotential and in the soft SUSY breaking 
potential. We generate output files compatible with CalcHep in order to use the public code 
micrOMEGAS_3.2 
with SuperIso (5^. 


52 . The R-physics observables are computed by interfacing the program 


The likelihood is constructed in a simple way. For measured quantities, we assume a Gaus¬ 
sian likelihood function with a variance given by combining in quadrature the theoretical and 
experimental variances. For observables for which only lower or upper limits are available, we 
use a likelihood modelled as a step function on the x% CL of the exclusion limit. The total 
likelihood function is then the product of the individual likelihoods associated to each experi¬ 
mental result. In order to save time in the sampling procedure, the slepton, chargino and gluino 
mass limits are, however, absorbed into the prior probability density functions: each parameter 
point generating a mass spectrum that violates one of these bounds is immediately discarded. 

Given the likelihood function, we sample the posterior probability density function with the 
MultiNest algorithm (54| - [^ . In order to cover all phenomenological interesting classes, we 
run separate chains that look either for hght(ish) EW-inos {jn^± < 900 GeV), light sleptons 
{nil < 600 GeV), or for light squark or gluinos {niq < 1.5 TeV or nig < 1.5 TeV). As for 
the choice of priors, we always take logarithmic priors on M^^niq, Aq^niHi while we use both 
logarithmic and flat priors for Mi, M 2 ,mL,niR,mN, Al, the sign of fi is fixed to -|-1 

(details on the prior ranges are provided in [^). In particular we perform two chains, one with 
log and one with flat priors, for each relevant data set: two chains for light EW-inos (these two 
data sets coincide with the ones used in j^), two chains for light sleptons and two chains for 
light squarks or gluinos. In each case, the other masses are left to vary freely from high to low 
values. The motivation for this is, as mentioned, to cover all potentially interesting cases; the 
results we will present in Section are for all chains combined together. 

The sampled points correspond to a 95% CL in volume of the posterior. (Since in this 
study we are not interested in statistical statements on the parameter space, we will however 
not exploit this feature.) The limits imposed by a step function are of course strictly obeyed by 
all scan points. Moreover, we have checked that none of the individual constraints implemented 
by a Gaussian gets a large pull in the final sample. In particular, BR(R —>■ A^q) and BR{Bs —> 
IjA iJi~) are in full agreement with the 95% CL experimental results for all points in the 

samples. 

Once the sampling of the parameter space according to the constraints in Table is com¬ 
pleted, all the points in the chains are confronted against the LHC Run 1 results by means of 
SModelS vl.0.1 as explained in the next subsection. 
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Figure 1: Illustration of “invisible compression” in SModelS vl.0.1. The decays of the heavier 
snentrino to the lighter one pins neutrinos are discarded in the hnal topology, leaving the Vr 2 
as an effective LSP. 


3.2 Deriving LHC constraints with SModelS 


SModelS vl.0.1 |34| - }36| is designed to decompose the signal of any arbitrary BSM spec¬ 
trum with a Z 2 symmetry into simplified model topologies and test it against the existing 


LHC bounds in the SMS context. SModelS vl.0.1 nses Pythia6.4 [^, NLL-fast |58| - [64| and 
PySLHA 1^, and includes a database of more than 60 SMS results from ATLAS and CMS. 
The decomposition procedure works “out of the box” for the MSSM-I-RN model with a sneu- 
trino LSP. Nonetheless some subtleties must be taken care of when processing the MSSM-I-RN 
scan points with SModelS vl.0.1. 

First, the inpnt to SModelS vl.0.1 can be simnlated events or an SLHA file contain¬ 
ing the fnll mass spectrum and decay tables as well as the SUSY prodnction cross sections, cr, 
in the format specified at [^. We choose the latter option. We use the MSSM-I-RN model 
implemented in micrOMEGAs 3.2 (see [^) to compute the decay branching ratios, B. 
The production cross sections for sleptons and sneutrinos [i.e. the sector modified with re¬ 
spect to the MSSM) are also compnted with micrOMEGAs 3.2. For all other prodnction 


processes, we nse the default SModelS vl.0.1 cross section calculator based on Pythia6.4 57 
and NLL-fast 58-64 . Electroweak cross sections are thus computed at leading order while 


strong productions are computed at NLO-I-NLL order. 

Given the information on a and B in the SLHA files, SModelS vl.0.1 computes a x B for 
each topology that occurs. Here, a topology is characterised by the SM particles originating 
from each vertex, and the mass vector of the SUSY particles in the decays. In order to avoid 
dealing with a large number of irrelevant processes, which is expensive in terms of computing 
time, topologies for which a x B < cTcut; with Ccut = 0.05 fb, are discarded. 

When dealing with an arbitrary spectrum of SUSY particles, it is possible that a part of the 
decay chain leads to completely invisible decays, e.g. a decay of a heavy sneutrino to a lighter 
one pins nentrinos in the current scenario. In snch cases, SModelS vl.0.1 compresses the 
invisible part of the decay chain as illnstrated in Fig. All decays to neutrinos appearing after 
the last visible decay are disregarded, yielding an “effective LSP” for the particular event, which 
can be different from the true LSP. This procedure is called “invisible compression”. Likewise, 
a neutralino may decay invisibly to a sneutrino and a neutrino; in this case the compressed 
topology resembles an MSSM topology. 

In addition, if the mass gap between mother and daughter particles is small, the decay 
products will be too soft to be detected at the LHC. This is taken care of by the so-called 


6 












“mass compression” in SModelS vl.0.1, discarding any SM particle that come from a vertex 
for which the mass splitting of the R-odd particles is less than a certain threshold. We use 5 
GeV as the minimum required mass difference for the decay products to be visible. 

Another comment is in order. The experimental constraints currently implemented in the 
SModelS vl.0.1 database require final states containing missing transverse energy (MET). 
This means that scenarios with long-lived particles (cr >10 mm) leading to signatures with 
displaced vertices or heavy charged particle tracks cannot be tested with SModelS vl.0.1. 
In the MSSM, this problem occurs, e.g., in wino-LSP scenarios where the yf is highly mass- 
degenerate with the and thus becomes long-lived. In the sneutrino LSP case, not only 
charginos can be long lived if the mass splitting with the sneutrino is small enough; other 
possibilities are, e.g., long-lived gluinos or stops, if they are the next-to-lightest SUSY particle 
(NLSP). We perform a detailed check of all input points to avoid the erroneous application of 
SMS limits to such cases. Points that have visible decays from long-lived particles or heavy 
charged particle tracks with cross sections larger than acut are discarded. (A brief discussion 
of such scenarios can be found in Appendix [B|) 

Once the decomposition into SMS topologies, including mass and invisible compression, is 
completed and the checks that the SMS results actually apply are passed, a given point is 
confronted against the SMS results in the SModelS vl.0.1 database. For each experimental 
constraint that exists, SModelS vl.0.1 reports among other things the analysis name, the 
Tx name identifying the topology]^ the predicted signal cross section for the point under con¬ 
sideration and the 95% CL experimental upper limit on it. Finally, the ratio r of the signal 
cross-section and the upper limit, r = cr (predicted) /cr (excluded), is given, where cr effectively 
means cr x 15 or the weight of the topology. A value of r > 1 means that the input model is 
likely excluded by the corresponding analysis. 

4 Results 

We now turn to analysing the impact of the LHC searches on the MSSM+RN parameter 
space. As explained in the previous Section, we here consider only points for which the SMS 
results apply, i. e. we discard points with non-prompt visible decays as well as points with long- 
lived charged particles. Scanning over the parameter space, we can then distinguish several 
cases: 

• the SMS results in principle apply but no SMS constraints actually exist for the specific 
topologies of the point — these points will be labelled as not tested^ 

• there exist (one or more) SMS results that test the specific topologies of the point but for 
each topology the total cr x H is below the corresponding 95% CL upper limit — these 
points will be considered as allowed; and 

• at least one topology has a a x B equal or above its 95% CL upper limit (r > 1) — these 
points will be considered as excluded. 

^The Tx names are explained in the SMS dictionary on http://smodels.hephy.at/wiki/SmsDictionary 
^This occurs if no simplified model result exists for the signal topologies of the point considered, but also if 
the mass vector of a topology lies outside that of the experimental constraint. Moreover, we include here also 
the points for which all signal topologies are discarded because of cr x .6 < Ucut- 
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Figure 2: For scan points that are excluded by the SMS limits, we show (in color) the breakdown 
of most constraining analyses in the vs. g mass plane. To illustrate the coverage of the 
parameter space, we also show (in grey) the not excluded or not tested points. 
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Figure 3: As Fig. [^but in the Vri vs. Xi mass plane. 


Let us start with the question which analyses are the most important ones for constraining 
the model. To this end, Figs. and show a breakdown of most constraining analyses in the 
versus g and i>ri versus Xi mass planes, respectively. Looking first at Fig.[^ we see that (the 
SMS interpretations of) the hadronic SUSY searches 


68 


■ 74 are constraining gluino masses up 
to about rrig ~ 1200 GeV and LSP masses up to about mc,^^ ~ 500 GeV. These searches mostly 
exclude points where either g —?• bb^, g —?• tt^ or g ^ qq^ decays are dominant, followed 


by an invisible decay of the neutralino, Xi Moreover, dilepton + MET searches 75,76 

exclude sneutrino LSP masses up to about 


210 GeV, independent of the gluino mass. 


The process that is constrained here is Drell-Yang production of xtXi followed by xf 
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Figure 4: Scatter plots of points for which SMS results apply. The top row shows the vs. 
g, the bottom row the vs. xf mass plane. In the panels on the left, the points excluded by 
the SMS constraints (red) are plotted on top of those which are not excluded (blue); in panels 
on the right this plotting order is inverted. Also shown (in grey) are the “not tested” points, 
for which no SMS constraints exist. 


(/ = e or /r), with the Pn + X decay being invisible (because of X being genuinely 

invisible or very soft). Consequently, in Fig. we see that chargino masses can be excluded 
up to about m~± ~ 440 GeV by the dilepton + MET limits. (There is also a small region of 
parameter space at low masses where T+r“ + MET 77 gives the strongest limit.) 

It is important to note here that the constraints on XiXi + MET actually stem 

from the —?• simplified model (and analogously for r+r" + MET), which has 

the opposite spin configuration than chargino-pair production followed by chargino decays into 
sneutrinos. The validity of applying the limits from the slepton searches to the case of chargino- 
pair production is discussed in Appendix |A] 

Also noteworthy is the fact that most of the excluded points in Figs. and have some 
grey points lying below them, which are not excluded or not tested at all. This is corroborated 
in Fig. 1^ where we present the summary of not tested, allowed and excluded points in the 
versus g and z>t-i versus xt mass planes. In the plots on the left, the excluded points (red) 
are plotted on on top of the allowed points (blue), while in the plots on the right this plotting 
order is inverted. As can be seen, only a small part of the parameter space can genuinely be 
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Figure 5: Allowed points in the vs. xf mass plane, with the color code indicating the 
wino/higgsino content of the Xi (|^7ii| = 1 means a pure wino while |J7ii| = 0 means a pure 
higgsino). 


excluded by the SMS results—over most of the regions where the SMS results are valid, there 
are almost always parameter combinations such that the limits can be avoided. 

For the dilepton signature originating from chargino-pair production, the chargino mixing 
plays an important role: wino-like charginos have a higher production cross section, and a 
higher branching fraction into Wn. The limits from + MET searches therefore mostly 
affect scenarios with wino-like xf^ while higgsino scenarios are much less constrained. For 
illustration see Fig. which shows the SMS-allowed points in the versus Xi mass plane— 
here the color map gives the size of the Uu entry of the chargino mixing matrix, indicating 
to the wino/higgsino content of the xt- As can be seen, in the region that is in principle 
constrained by the SMS results the surviving points feature yfs that have a large higgsino 
admixture (|17ii| ^ 0.5). These points have a lower xtxi production cross section and the xt 
decays preferably into Ti>r^ since the higgsino decay to e, /r is Yukawa suppressed; r+r" -l-MET 
is however a more difficult signature experimentally and thus only constrains a small strip at 
low P mass, cf. the purple points in Eig. 

Missing topologies 

The next question to ask is which are the most important signatures not covered by SMS 
results. Such information can be used to improve on the interpretation of the LHC searches for 
new physics. We call these uncovered signatures “missing topologies”. For any point passed 
through SModelS vl.0.1, we keep up to ten missing topologies sorted by their a x B. To avoid 
double counting, missing topologies are evaluated after mass and invisible compressions. The 
total weight is computed by summing over all diagrams giving the same topology, i.e. ignoring 
the mass vector of the SUSY states involved. Moreover, I = e, fj, lepton flavors appearing in 
the final state are summed over (light quark flavors are always summed over). In the following, 
we only consider MSSM-I-RN scan points which are not excluded, and we demand that missing 
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Figure 6: Missing topologies with highest a x B in the Un vs. g mass plane. 


topologies have a x B > 1 fb. The results can be presented in two ways, either by showing 
the most frequent missing topologies in a certain parameter space, or by selecting for each 
parameter point the missing topology with the highest cross section. 

We choose the latter approach to show in Fig. the missing topologies in the sneutrino- vs. 
gluino-mass plane. The various processes are denoted in the bracket notation of SModelS, ex¬ 
plained in 34 . The structure is [branchl, branch2] for the decay chains (“branches”) of the 


two initially produced SUSY particles; each branch contains inner brackets for each vertex, con¬ 
taining in turn the lists of outgoing standard model particles. Thus, [ [ [b, b] ] , [ [ j et , j et] ] ] 


denotes gluino-pair production with one gluino decaying into bb + MET (via g —?• 66%°, 
uu) and the other one into qq + MET. Likewise, [[[jet, jet]] , [[jet, jet] , [1]]] 
denotes gluino-pair production with the first gluino decaying via g qqx^, and the 

other one via g qq'x^, 

It is apparent that many points with gluino masses below about 1.2 TeV, for which the 
LHC searches should have good sensitivity, are not excluded by the SMS results because they 
feature “mixed topologies”, where the two pair-produced gluinos undergo different decays {e.g. 
one gluino decaying into bb and the other one into light jets). Since the SMS results for pair- 
produced sparticles always assume two identical branches, these cases cannot be constrained by 
SModelS vl. 0.1. Moreover, hadronic final states with additional leptons, as they arise from 
gluino decays into charginos and the chargino decaying further into a charged lepton (e, /r or 
r) plus the LSP, do not have any SMS equivalent. Einally, there are no SMS results available 
for g tbxf, no matter of whether the chargino has any visible decays. 

It is also worth noting that over a large part of the parameter space single lepton -|- MET 
([[],[[!]]] in bracket notation) is the most important missing topology. This signature arises 


®More generically, [[[b,b]] , [[jet,jet] ,[!]]] denotes production of XY with X undergoing a 1-step 
decay chain, X ^ bh -\- MET (brcinchl= [vertexl] = [[b,b] ] ) and Y undergoing a 2-step decay chain, Y —> 
qq + Z qq +1 + MET (brainch2= [vertexl,vertex2] = [ [jet, jet] , [1] ]); X can be different from Y or both 
can be the same. 
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Figure 7: Missing topologies with highest a x B va the vs. Xi mass plane. 


from Xi production; its importance is corroborated in Fig. where one can see that it is 
indeed dominating the whole sneutrino- vs. chargino-mass plane. (There are also cases where 
single W + MET is dominant.) The cross section for single lepton + MET production, shown 
in Fig. can be very large and should give important additional constraints on the model. 
While searches for single lepton + MET were performed by both ATLAS and CMS (7^ , 
unfortunately no suitable SMS interpretation exists for these analyses. It would be extremely 
interesting if the experimental collaborations provided upper limit maps and/or efficiency maps 
for their single lepton + MET analyses in the context of a chargino-sneutrino simplified model. 

Having both light EW-inos and light staus can generate decay chains with more ‘exotic’ 
signatures, in particular xf^j followed by xf and x^ —>■ —>■ t^W^p- 


'Tl ■ 


This appears as [ [ [nu] , [W] ] , [ [ta] , [W] ] ] (yellow points) in Eig. [^and is interesting because 
the x° decay produces with the same rate t~^W~ and together with the chargino decay 

this gives rise to a same-sign W signature, + MET. 

Before proceeding it is instructive to take another look at the missing topologies arising 
from EW-ino and slepton production, but this time ordered by their frequency of occurrence. 
This is done in Eig. Not surprisingly we see that besides single lepton (e or ;u), single r is an 
important signature. Although it is less clean experimentally, the relative weight of single e, /r or 
T -h MET might potentially give information on the mass pattern of the mostly RH sneutrinos. 
Another important class of “missing topologies” are different-flavor dileptons ([[[!]] , [ [ta] ] ] 
and [[[!]] , [[!’]]]). Different-flavor dileptons -|- MET have in principle been considered by 
ATLAS and CMS in the context of chargino-pair production in the MSSM with the charginos 
decaying either into kF^*^X° 75 or into/z^x)* via on-shell sleptons/sneutrinos 75,76 . However, 


the associated SMS limits do not apply to the sneutrino LSP case for various reasons. For 
example, the leptons from xf ^ are generally softer than those from xf —?• l^Pi decays 

(for the same xf and LSP masses) because of the additional neutrinos in the W decay. The 
limits for the XiXi 2 x Ip{ov pI) —>■ 2 x IpXi simplified model are also not applicable because 
they involve an additional intermediate mass scale. 
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Figure 8: Cross sections a x for the single lepton + MET missing topology for not excluded 
or not tested points. 
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Figure 9: Missing topologies with a xB > 1 fb in the sneutrino- vs. chargino-mass plane ordered 
by frequency of occurrence. The ordering is from top to bottom in the legend, with single tau 
being the most frequent missing topology, followed by single lepton (/ = e,p), lepton-tau, and 
so on. “Other topologies” are shown on top of the legend without considering their total count 
(however, each single one of them is less frequent than any of the topologies denoted explicitly). 


Finally, the [ [ [W] ] , [ [W] , [ta] ] ] topology again gives rise to same-sign IF’s, see the red tri¬ 
angles in Fig.|^ Similarly it is possible to have same sign r’s arising from [ [ [W] , [ta] ] , [ [ta] ] ] 
(black stars). In this case, after x^xf production, the decay chain is Xi 
and xf 
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Most constraining SMS results 
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Figure 10: Complementarity of LHC and direct DM detection experiments. The panel on the 
left shows SMS allowed, excluded and not tested points in the plane of vs. . The panel 
on the right shows the breakdown of most constraining analyses for the points that are excluded 
by the SMS limits (for the sake of comparison, the allowed points are shown in grey). In both 
panels, the solid magenta lines and the dashed blue lines are the current exclusion limit by 
LUX and the forecasted sensitivity of XENONIT experiment respectively, while the dashed 
light green line corresponds to the predicted neutrino coherent scattering on nuclei. 


Complementarity with direct DM searches 


Let us finally turn to the complementarity of LHC and direct DM searches—recall that all 
points in our scans are consistent with DM constraints, as described in Table In Fig. left 
panel, we plot the allowed (gray), excluded (red) and not tested points (cyan) as a function 
of the sneutrino mass and the SI scattering cross section. In the same plot we also show the 
forecasted sensitivity of XENONIT after two years of scientific run and the predicted value 
for neutrino coherent scattering on nuclei , which can be an irreducible background for direct 
detection experiments. From this plot, the complementarity between the two type of searches 
is striking. Points with a SI elastic cross section well below the neutrino background, and hence 
not detectable by direct detection experiments, are already excluded by SMS results. On the 
other hand, a bulk of points allowed (or even more interestingly, not tested) by SMS results is 
well in the reach of XENONIT, expected to start running in 2015. Notice however that there 
still exist combinations of parameters that allow sneutrino DM to escape both direct detection 
and LHC searches, represented by the cyan points below the neutrino background curve. In the 
MSSM+RN, DM direct searches are basically sensitive to the mass of the LSP and its couplings 
with the Higgs and Z bosons. The rest of the SUSY mass spectrum is not relevant. This is 
different with respect to the MSSM with the neutralino LSP, where the interaction with the 
quarks is mediated as well by squarks on t-channel. This is clearly visible in the right panel of 
Fig. which shows the most constraining SMS analyses. In Figs. and these SMS analyses 
are typically correlated with the gluino or chargino mass, while now they are scattered all over 
plane. 


the versus m,-. 
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Figure 11: On the left allowed (grey), excluded (red) and not tested (blue and cyan) points are 
shown in the plane of sneutrino mass versus mixing angle. The subset of points with exceedingly 
small a X B at the LHC but in the reach of XENONIT is visualised in blue. On the right, we 
show the SMS allowed points in the sneutrino mass versus mixing angle plane, subdivided in 
blue points, which are in the reach of XENONIT and in light green (light grey) points, which 
are above (below) the neutrino background. 


The same set of allowed, excluded and not tested points are plotted as a function of the 
sneutrino mixing angle in the left panel of Eig. [T^ The bulk of not tested points in the reach of 
XENONIT (dark blue points) has, as expected, relatively large mixing angles, corresponding to 
sizeable contributions from Z boson exchange to the SI scattering cross section. Excluded red 
points are scattered everywhere in the sin dp vs. plane and probe also very RH sneutrinos. 
In the right panel of Fig. [m we see that among the allowed points, XENONIT can constrain 
a large portion of the sneutrino parameter space, while the very RH sneutrinos will remain 
inaccessible to future direct detection detectors. In general the points with negligible mixing 
angles have as LSP and the neutralino as NLSP, which tends to be almost degenerate with 
chargino. The relic density is then actually achieved by co-annihilation of neutralino-chargino 
and then communicated to the mostly sterile LSP (see for details). Such scenarios are very 

difficult to test. 
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5 Conclusions 


Scenarios with a sneutrino as the LSP are an interesting alternative to MSSM models with 
neutralino LSPs. Indeed in SUSY models with a RH neutrino superfield (MSSM+RN) the 
fermionic field contributes to neutrino masses while the scalar field contributes to the DM 
candidate, which is a mixed, however mostly RH, sneutrino. 

The collider phenomenology of the MSSM+RN can be quite different from the typical MSSM 
case. It is therefore interesting and relevant to ask how the SUSY search results from Run 1 of 
the LHC, which were mostly designed with the MSSM in mind, constrain sneutrino LSP sce¬ 
narios. To address this question, we used SModelS vl.0.1 for testing the MSSM+RN against 
more than 60 results from CMS and ATLAS searches in the context of so-called Simplified 
Model Spectra (SMS). More precisely, by considering the model parameter space where the 
sneutrino is a good DM candidate compatible with all current constraints, we assessed 1.) the 
constraining power of the current SMS results on such scenarios and 2.) the most relevant 
signatures not covered by the SMS approach. 

Concerning point 1.), we found that the dilepton + MET searches are among the most rele¬ 
vant ones, constraining sneutrino masses up to about 210 GeV and mostly wino-like charginos 
up to m~± ~ 440 GeV. It is important to note here that this amounts to re-interpreting the 

ATLAS and CMS searches for pp —?• in terms of pp —?• XiXi (the 

validity of this is discussed in Appendix]^. Hadronic SUSY searches exclude gluinos masses up 
to rrig ~ 1200 GeV and LSP masses up to rric,^ ~ 500 GeV. Nonetheless in general we find that 
only a very limited portion of the parameter space can be properly excluded by SMS results. 
For most points in the {mg,mp^) or (m~±,mj>J planes there exist parameter combinations that 
allow to avoid all limits. Indeed, most of the parameter space is either allowed (SMS constraints 
exist for the specific topologies of the point but all a x H of these topologies are below their 
95% CL upper limits) or not tested at all (there are no existing SMS constraints for the specific 
topologies of the point or each topology has a, a x B which is smaller than 1 event at LHC 
Run 1). Direct DM searches are complementary to the SMS constraints: many points that 
are not tested by SMS results can potentially be excluded by XENONIT. Vice versa, points 
well below the neutrino background, hence not reachable by future DM detectors, are already 
excluded by SMS results. 

The second main result of this paper concerns point 2.), i.e. the study of the allowed points 
in terms of missing topologies. In the hadronic sector, pair-produced gluinos with masses well in 
the reach of LHC Run 1 are not constrained because they feature one or more of the following: 

• additional leptons: since the gluino cannot directly decay into the sneutrino LSP, the 
hadronic final state is often accompanied by leptons; 

• mixed topologies: each of the pair-produced gluinos undergoes a different decay; 

• the gluinos decay into tb final states. 

None of these possibilities are covered by the current SMS results. Note here that the last 
two items are also common in the MSSM, as described in [^. For EW production, missing 
topologies include: 

• single leptons; 
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• single lys; 

• different-flavour opposite-sign leptons; 

• same-sign Ws or same-sign tans (accompanied respectively by additional leptons/taus, 
or ffls). 

While such signatures have been searched for by the SUSY and/or exotics groups in ATLAS 
and CMS, the results do not exist in terms of appropriate SMS interpretations. Such an SMS 
interpretation would be very interesting in particular for the mono-lepton -|- MET case, which 
promises to have a considerable impact for constraining the MSSM-I-RN modelj^ 

A final comment is in order. While the SMS approach is very convenient for the charac¬ 
terisation of new physics signatures and vast surveys of parameter spaces, it clearly has its 
limitations. Given the high interest in non-standard SUSY (and other new physics) scenar¬ 
ios, we urge the experimental collaborations to document their analyses in a way that they 
can conveniently be re-casted in public simulation frameworks like CheckMATE or the 

. (See also the recommendations in and 
This would allow to go beyond the limitations of SMS approach and give a much more rigorous 
assessment of the constraints in a large variety of new physics models, including the sneutrino 
DM scenario discussed in this paper. Unfortunately we are still a long way from this. 


84 in this context). 
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®This could be done analogous to the existing xt^ ixt W^XiJ X 2 ^^Xi) simplified models that are 
already assessed by the ATLAS and CMS SUSY groups, but with the chargino decaying to 100% into 
and the neutralino decaying 100% into lyiyi. However, since the chargino and neutralino masses need not be 
degenerate, we propose to consider as a first step xtxi production followed by xf ^ Xi 

The cross section upper limits should be provided in the chargino- versus sneutrino mass plane for different 
neutralino masses, for the cases I = e,iJ, and I = r, and if computationally feasible also for I = e,ii,T assuming 
equal rates. 
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A Validity of slepton search results for chargino-pair 
production with decay into lepton+sneutrino 

In the spirit of the SModelS vl.0.1 philosophy, we apply SMS constraints for the + 
MET topology, obtained in the context of pair production of charged sleptons, pp 
followed by —>■ l^^Xi fo the case of chargino-pair production pp —>• XiXi followed by 
xf despite the opposite spin configuration. This can only be valid if the signal 

selection efficiencies in both scenarios are comparabfe. To test this assumption, we use the 
recast code 85 for ATLAS search in final states with two leptons and missing transverse 


momentum, ATLAS-SUSY-2013-11 [^, which is available in the framework of the Mad- 
AnalysisS “Public Analysis Database” |^. We consider two benchmark scenarios in the 
simplified-model spirit, an MSSM one with {mj±, m^o) = (270, 100) GeV and an MSSM-I-RN 


m~±, mc ,^) = (270, 100) GeV. Events are generated with MadGraph 5 


one with 

Pythia 6.4 |57| and then passed through Delphes 3 


86 


87 


and 


for the simulation of the detector 
effects]^ Eor simplicity, in the following we restrict our study to pair-production of selectrons 
for the MSSM case, and pair-production of charginos decaying exclusively via electrons in the 
MSSM-I-RN case. 

The event selection requires two opposite sign (OS), same flavor (SE) leptons with high 
transverse momentum, concretely pr > 35 GeV and pr > 20 GeVp] Eigure [T^ compares the 
Pt distributions in the two benchmark scenarios, in the left panel for the harder electron, ei, 
in the right panel for the second electron, 62 . The bin sizes are chosen such that the first bin 
corresponds to the events that do not pass the pt > 35 GeV (left panel) or pt > 20 GeV (right 
panel) requirement. We see that the electrons originating from selectron-pair production tend 
to be harder than those originating from chargino-pair production. 

The analysis further requires the invariant mass of the lepton pair to be outside the Z 
window, and rs and jets are vetoed. Einally, three signal regions are defined by thresholds 
on the mT 2 (“stransverse mass”) variable that is used for reducing the tt and Wt 

backgrounds: mT 2 > 90, > 120 and > 150 GeV. The mT 2 distributions after the preselection 
cuts are shown in Eig. [T^ It can be seen that the distributions intersect around the minimum 
required value of 771^2 = 90 GeV; events with electrons originating from chargino decays are 
more likely to pass this cut. 

To see the net effect on the signal efficiencies, Tabie[^ shows the complete cut-flow compar¬ 
ison for the two benchmark scenarios. As expected, differences arise in the first cut, selecting 
high Pt OS lepton pairs, and when applying the lower bounds for mr 2 - Because of the softer 
Pt distribution in case of chargino production-|-decay, there are fewer events passing the first 
cut for this scenario. However, the opposite is true for the mT 2 cut. Ultimatefy, the efficiencies 
are comparable in all signal regions, and even somewhat higher for the MSSM-I-RN scenario. 

To check that this is still true closer to the kinematic edge, we reproduce the cut-flows for 
a second set of benchmark scenarios with an LSP mass of 200 GeV. As can be seen in Table 
we And a similar behaviour in this case. We conclude that we can safely apply the SMS upper 


^Note that for the reconstruction of events with a sneutrino LSP it is necessary to define the sneutrino as 
MET, by adding a corresponding EnergyFraction entry in the Delphes card. 

®We consider here only the part of the analysis that is relevant for the SMS result used to constrain the 
sneutrino LSP scenario in Section 4. 
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Figure 12: Comparison of the pt distributions of electrons originating from selectron decays 
in the MSSM and from chargino decays in MSSM+RN, at the level of reconstructed events. 
The benchmark scenarios used are (m;-±, m^oj = (270, 100) GeV for the MSSM case and 
(m~±, = (270, 100) GeV for the MSSM+RN case. See text for details. 



Figure 13: Comparison of the mT 2 distributions for the two benchmark scenarios after all 
preselection cuts. 

limits given by the experimental collaborations in the context of slepton-pair production in the 
MSSM to constrain chargino-pair production followed by decays into lui in the MSSM+RN. 
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Table 2: Comparison of the cut-flows for pp ^ ee —>■ e+e XiXi pp ^ Xi^Xi “t e+e vivi 
with (m[±, m^o) = (270, 100) GeV and mi,^) = (270, 100) GeV, respectively. 


Cut 

Slepton production 

Ghargino production 

Common preselection 

Initial number of events 

50000 

50000 

2 OS leptons 

35133 

33464 

mu > 20 GeV 

35038 

33337 

r veto 

35007 

33318 

ee leptons 

35007 

33318 

jet veto 

20176 

19942 

Z veto 

19380 

18984 

Different mT 2 regions 

mT 2 > 90 GeV 

11346 

11594 

771x2 >120 GeV 

8520 

8828 

771x2 >150 GeV 

5723 

5926 


Table 3: As Tablebut for {mi±, 'm^o) = (270, 200) GeV and = (270, 200) GeV. 


Cut 

Slepton production 

Ghargino production 

Common preselection 

Initial number of events 

50000 

50000 

2 OS leptons 

29291 

27244 

mu > 20 GeV 

29082 

26964 

T veto 

29050 

26956 

ee leptons 

29050 

26956 

jet veto 

16834 

16114 

Z veto 

15281 

14025 

Different 771^2 regions 

771x2 > 90 GeV 

3028 

3198 

771x2 > 120 GeV 

85 

140 

771x2 >150 GeV 

0 

0 
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B Lifetimes of long-lived particles 


As mentioned in Sec. |3.2[ a considerable number of the scan points comprise long-lived spar- 
ticles. These occur mostly when enforcing light gluinos or squarks; in this case about 30 % of 
the points feature long-lived particles, while the fraction is below 1 % without this constraint. 
The long-lived particles are predominantly gluinos (85 %), mostly in the case where it is the 
NLSP, and in a few points where Xi is slightly (up to about 50 GeV) lighter than the gluino. 
Apart from that we find points with long-lived stops or staus in case they are the NLSP, as 
well as single points with long-lived charginos. Here we will focus on the long-lived gluinos and 


stops, long-lived staus have been discussed before in 27 



m;j [GeV] [GeV] 


Figure 14: Lifetimes cr in [m] for long-lived gluinos, the color code indicates the LSP mass 
(left) and the sneutrino mixing angle (right). 

In the MSSM long-lived gluinos appear when all squarks are extremely heavy, e.g. in split- 
SUSY scenarios. In case of the MSSM-I-RN with a sneutrino LSP additional causes come into 
play. If the gluino is the NLSP, its decay will proceed only via virtual squarks and gauginos, 
yielding an effective four body decay, g —>■ gqi'P (virtual q and x°) or g ^ qq'li) (virtual q and 
X^). The gluino lifetime will therefore depend not only on the squark mass, but also on the 
gaugino masses and mixings, as well as the sneutrino mixing angle. Meta-stable gluinos can 
thus appear even if the squarks are not completely decoupled. The gluino lifetime as a function 
of its mass is shown in Fig. The left plot illustrates the depencence on the sneutrino mass, 
the right plot the dependence on the sneutrino mixing. We can distinguish two general regions. 
First, we observe an exponential dependence of the lifetime on the gluino mass for decay lengths 
of 10 mm up to 10^ m. Here the lifetime is largely independent of the sneutrino mass. Moreover 
lifetimes at constant gluino masses are longer for heavier squarks and gauginos. In this region 
we generally find large mixing angles sin dp, but heavy gauginos and squarks. Points with very 
small mixing angles may also appear in this region, in the case that the mass of the lightest 
neutralino is below the gluino mass. The second region, with lifetimes longer than 10^ m, and 
up to 10^^ m, shows a very different behaviour. We can see a clear correlation between gluino 
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and sneutrino masses in this region, with longer lifetimes found for smaller mass splittings. 
The lifetimes moreover increase when going to very small sneutrino mixing angles, with the 
maximum lifetimes achieved for sin going to zero. 

Likewise, if the stop is the NLSF|^ and has a small mass difference with the sneutrino, it 
can be long-lived, see Fig. [T^ As seen for the gluinos, the lifetime depends strongly on the 
sneutrino mixing. 



Figure 15: Lifetimes cr in [m] for long-lived stops, the color code indicates the LSP mass (left) 
and the sneutrino mixing angle (right). 


Both long lived gluinos and long lived stops can be constrained by searches for R-hadrons, 
see |9l||^ for R-hadrons escaping the detector, for stopped R-hadrons, or 94 for metastable 
gluinos decaying in flight inside the detector. However, large uncertainties arise from modeling 
both the hadronisation and the strong interaction of the R-hadron with the detector. Therefore 
no collider constraints on long-lived sparticles have been included. 

Additionally, cosmological constraints become important for gluino lifetimes of about 100 s 
(10^° m) j^. Lifetimes of that order would affect the fraction of heavy nuclei produced dur¬ 
ing the Big Bang nucleosynthesis. Longer lifetimes can further be constrained by searches for 
diffuse gamma ray background, distortions in the CMBR and heavy isotopes. 


®If the stop mass is close to the gluino mass, both stop and gluino may be long-lived. 
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